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dimensions. In this framework we observe an area law for spatial Wilson loops in four 
and five dimensional supersymmetric Yang-Mills at finite temperature. This can be 
interpreted as the area law of ordinary Wilson loops in three and four dimensional 
non-supersymmetric gauge theories at zero temperature which indicates confinement 
in these theories. Furthermore, we show that super Yang Mills theories with 16 
supersymmetries at finite temperature do not admit phase transitions between the 
weakly coupled super Yang Mills and supergravity regimes. This result is derived 
by analyzing the entropy and specific heat of those systems as well as by computing 
ordinary Wilson loops at finite temperature. The calculation of the entropy was 
carried out in all different regimes and indicates that there is no first order phase 
transition in these systems. For the same theories at zero temperature we also 
compute the dependence of the quark anti-quark potential on the separating distance. 
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1. Introduction 

In the last couple of months it has become clear that supergravity is a useful tool 
to study the large N limit of field theories [1]. Related work and many new results 
can be found in refs. [2] - [35]. In the present paper we follow this approach and 
consider phase transitions, Wilson loops and confinement in super symmetric and non- 
super symmetric field theories using supergravity. We use the corresponding extremal 
and near extremal supergravity solutions in the decoupling limit [1, 10] and study 
Wilson loops. To consider the supersymmetric theories at zero (finite) temperature 
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we follow [37, 38] ([39, 40]) and use the Nambu-Goto action in the background of an 
extremal (near-extremal) D-brane solution to calculate the space-time Wilson loop. 
We find that for the field theories with maximal supersymmetries there is no finite 
temperature phase transition between the SYM and SUGRA regimes [10]. This we 
also check via a direct entropy consideration which shows that the entropy matches 
(up to numerical coefficients which we do not calculate) across the different domains. 
Since the entropy seems to match in all the different regimes of the theory it is most 
likely that there is no first order phase transition. This does not exclude higher 
order phase transitions which are known to occur in various cases as we enter the 
eleven dimensional M theory regime. In those cases a "localization" phase transition 
takes place as discussed in [10]. To consider non-supersymmetric theories at zero 
temperature we follow [46] and consider the spatial Wilson loop in the background 
of Euclidean near-extremal Dp-brane solutions. When the spatial size is much larger 
then 1/T (where T is the Hawking temperature of the near-extremal solution) the 
effective low energy theory reduces effectively to a p dimensional non-supersymmetric 
theory. Therefore, the spatial Wilson loop gives us the energy between a quark and an 
anti-quark of the p dimensional theory at zero temperature. Using this approach we 
find an area law behavior for non-supersymmetric YM in three and four dimensions 
in the large N limit. 



The paper is organized as follows: In section 2 we briefly describe following [10] 
the supergravity solutions of N coincident Dp-branes in the field theory limit that 
was introduced in [1]. In Section 3 we use the supergravity approach to confirm the 
confinement behavior of non-supersymmetric YM theory in three (four) dimensions. 
We consider the Euclidean theory of N D3- (D4-) branes with ,R 3 x S 1 (i? 4 x S 1 ) 
world volume. In the limit of small radius of the S 1 circle, imposing appropriate 
boundary conditions, we end up with a three (four) dimensional pure YM theory. 
In this setup we compute the spatial Wilson loop and show that it admits an area 
law behavior. Section 4 is devoted to analyzing the entropy and the specific heat 
of the various theories with 16 supersymmetries. It is shown that there is no phase 
transition in the SYM, supergravity and other domains of the Dp-brane systems. In 
section 5 we present the derivation of ordinary Wilson lines (along one space and one 
time directions) for these theories at zero temperature. This includes the derivation 
in the ten dimensional supergravity regime, as well as in the extensions to other 
energy domains. In section 6 the finite temperature behavior of these systems is 
deduced from the supergravity description. It is shown that for p = 1 ... 4 they 
admit a similar behavior as the one discovered in [39], [40]. 
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2. A very brief review of the theories with 16 supersymme- 
tries 

In [10] systems of N coincident extremal Dp-branes where analyzed in the decoupling 
limit [1] 

g \ = (2n) p - 2 g s a' {p ~ 3)/2 = fixed, a -> 0, U = ^ = fixed, (2.1) 

a 

where g s = e^°° , and g YM is the coupling constant of the p + 1 dimensional U (N) 
SYM theory ( with sixteen supercharges) that lives on the N Dp-branes. In the SYM 
picture U corresponds to finite Higgs expectation value associated with a U (N+l) — > 
U(N) x U(l) symmetry breaking. The effective coupling of the SYM theory is 
Qeff = 9ymNU p ~ 3 . Perturbation theory can be trusted in the region 

9eff « 1. (2.2) 

The type II supergravity solution describing N coincident extremal Dp-branes can 
be trusted if the curvature in string units and the effective string coupling are small. 
These conditions yield 

1 < g 2 eff < . (2.3) 

which translate for p < 3 to the following range of the energy scale U 

(^Ar)V(3-p) 7V 4/(p-3)(7-p) <<; u <<; {g 2 YMN) l/(3- P ) {2A) 

For p > 3 the signs are replaced with ^> ones. In the supergravity description U 
is the radial coordinate. 

In general, there are other regions (which will play a role below) which take over 
when the dilaton becomes large. We do not describe these regions here since they 
depend on p and cannot be describes in this concise manner. For details see [10]. 

3. Area law 

In this section we use the supergravity approach to field theory and find confinement 
in non-supersymmetric theories in three and four dimensions. 

3.1 Area law for three dimensional YM theory 

We consider the Wilson loop along two space directions in the case of the near 
extremal D3 brane solution. We will find that it shows area law behavior as in 
the R 3 case in ref. [46]. We shall take one direction, Y, to be large and the other 
direction, L, to be finite. In the limit Y — > oo we consider configurations which are 
invariant under translation in the Y direction. 
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To describe the theory at finite temperature we go to a Euclidean description and 
compactify the time direction on a circle with period (3 = T _1 . For large temperature 
this circle is small and the theory becomes effectively the Euclidean description of a 
three dimensional field theory with gauge coupling given by 

9ym 3 = 9ym 4 T ■ (3.1) 

Since we choose boundary conditions on the circle such that supersymmetry is bro- 
ken, the fermions and scalar fields are heavy with masses of the order T and gy M4 T, 
respectively [46]. Therefore, at large distances, L ^> f3, we obtain zero temperature 
non-supersymmetric Yang-Mills theory in three dimensions. Confinement and area 
law for the Wilson loop are expected. We will derive this momentarily by calculating 
a spatial Wilson loop of the "compactified" four dimensional theory. Since one of 
the space directions becomes the Euclidean time direction of the three dimensional 
theory this Wilson line is an order parameter of the theory. Note that the area law 
behavior, found in this section, does not imply confinement in the 3+1 dimensional 
theory with temperature. 

Confinement is expected to appear in the limit where LT ^> 1 i.e. the distance 
between the quarks is much larger than the size of the circle. On the other hand if 
we consider small distances we are back in the zero temperature 3+1 dimensional 
theory studied in [37]. Thus, our description interpolates between confinement in 
2+1 dimensions (LT 3> 1) and Coulomb behavior in M = 4 YM in 3+1 dimensions 
(LT < 1). 

The metric of near extremal D3 branes in the large N limit is 



ds 2 = a' { ^[-f(U)dt 2 + dx 2 } + R 2 f(Ur^ + R'dQl 

f(U) = 1 - U*/U 4 (3.2) 

R 2 = ^g~N, U* = yTrV/x , 

where /x is the energy density. Thus, the relevant action for the spatial Wilson loop 
is 

„ Y 



jdx, 



~ 2tt J \\ i? 4 1 - U^/U 
The distance between the quark and the anti-quark is 

R 2 f°° dy 



u A (d x uy . . 

+ n 77I777I ' (3-3) 



J u Ji 



(3.4) 

V-l)(y 4 -A) 

where A = U^/Uq < 1 and Uq is the minimal value of U. Notice that in the limit 



Uq — > Ut (A — > 1) we get L 




1|+^^, (3.5) 

7T 
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where, as was explained in [39], the subtraction is at the horizon. We are after the 
large L limit behavior. Thus we need to take the limit A — > 1. In this limit the main 
contribution to the integrals in (3.4) and (3.5) comes from the region near y — 1. 
Therefore, we get for large L 

E = T QCD L. (3.6) 

The tension of the QCD string is 

Tqcd = \R 2 T\ (3.7) 

where we have used the relation [39] Ut = nR 2 T. 

It is important to emphasize that the theory for which we obtain the area law 
is not YM in three dimensions but Af = 4 at four dimensions compactified on a 
circle. The string tension which we derive, (3.7), "knows" about the four dimensional 
origin of the theory. The reason is as follows. From (3.7) we see that the mass of 
the excitations of the QCD strings is M s = RT. The mass associated with the 
compactification is M c = T. To trust the supergravity solution we need R >> 1 [1]. 
Thus M s ^> M c and so the QCD strings probe distances which are much smaller 
then the radius of compactification. 

Some speculations about string theory at large curvature can be made. For 
YM in three dimensions we expect to have Tqcd ~ 9% which in our notation is 
~ R 8 T 2 and not ~ R 2 T 2 as in (3.7). The system which we are considering here 
is a good approximation to YM in three dimensions if R <C 1 (because then the 
masses of the excitations of the QCD string is lighter then the masses associated with 
the compactification). The supergravity solution cannot be trusted in that region. 
However, since the solution is near-extremal it makes some sense to speculate that 
the exact string theory solution has the same form but with corrected harmonic 
functions. In that case in the limit which we are studying we will get for the exact 
string theory solution an AdS space times a sphere where the radius of the AdS is 
not R 2 but a function of R 2 . The fact that for YM in 3d we expect to get R s for 
R 2 <C 1 implies that the exact string theory solution interpolates between R 2 at large 
R 2 and R s at small R 2 . 

3.2 Area law in four dimensional YM theory 

The approach of the previous section to confinement can be generalized to obtain 
confinement in four dimensions from supergravity. We need to consider the super- 
gravity solution of near-extremal D4-brane in the decoupling limit. A D4-brane is 
described in M-theory as a wrapped M5-brane so from the point of view of M-theory 
we relate the near-extremal solution of M5-brane to confinement in four dimensions 
as was suggested in sec. 4 of [46]. 
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The near-extremal solution of D4-branes in the decoupling limit is [10] 



a 



R 3/2 

-K4 



3/2 



1 - U T /U 3 )dx 2 + dx\ + ... + cfa^) + 



dU 2 + R 3 d /2 Vudn 2 



(3.8) 



(27T) 



i? 3 



3\ 1/4 



where i? 3 ^" 



■4 — 9§y 7^ an( i #5 i s the 5D SYM coupling constant. 



We would like to study the spatial Wilson loop in the region where LT ^> 1. In 
this region the effective description is via a non-sup ersymmetric YM theory in four 
dimensions with coupling constant 



9ym = fir. 



(3.9) 



Unlike the supergravity solution which was used in the previous section the super- 
gravity solution (3.8), which we use here, cannot be trusted for arbitrary U. L, is 
related to U by L ~ 1/U and hence it is also bounded. 

Before we perform the calculation of the spatial Wilson line let us first find the 
upper bound on LT and the bounds on g YM and g\^ = g YM N. The restrictions 
on U and hence on Ut, are such that the curvature in string units and the effective 
string coupling are small. The result of these restrictions is (2.4) [10] 



1 TT N 1 ' 3 

-r^r-y < U T « 5-. 

Ngl gl 
Therefore, the supergravity solution can be trusted only for distances 



Ngl > L > 



(3.10) 



(3.11) 



To find the bound on T we use the relation between T and Ut (the temperature can 
be obtained from (3.8) andT = ^ 



T = 



U=U T 

3 



2V^Ng, 



U T , 



to get 



1 nn 1 

< T < 



Ng 2 



N^gl 



(3.12) 



(3.13) 



From (3.9) we find that the four dimensional coupling constant is bounded by 

Jf « 9ym « ^73 • (3-14) 
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We see, therefore, that in the large N limit qym must go to zero. For the four 
dimensional effective coupling, g e ff, we have 

1 < g 2 en < iV 2 / 3 . (3.15) 

Thus the effective four dimensional coupling constant g e ff must be large otherwise 
we cannot trust the supergravity description. Finally we turn to the bound for TL. 
To be able to use the supergravity results described below we need to find a region 
where TL > 1. From (3.11) and (3.13) we get 

iV 2 / 3 » TL » j^. (3.16) 

Therefore, in the large N limit there is a region for which we can trust our results. 
Note that unlike the 3d case, considered in the previous section, in the 4d case, for 
any finite N, L is bounded. 

Let us now derive the area law behavior. The action for the string in this case is 



— f dx* 
2tt J 



IP d x U 2 



Using the same manipulations as in [37] we get 

dy 




£ = -/ dy\ , * -i + gZ^g°, (3.18) 

D!// :i A;) J 

jj3 

where A 4 = For TL > 1 we have (Ut — U )/U <C 1 and the integrals are 

dominated by the region close to y — 1. Therefore, as in the previous section, we get 

E = T YM L (3.19) 

where the string tension is 1 

T YM = ^g YM NT 2 . (3.20) 

This agrees with known large iV results if T is identified with Aqcd up to a N 
independent constant factor. Again, as in the previous section, in the region where 
we can trust the supergravity solution, the QCD string can probe distances which 
are much smaller than the compactification radius and hence it "knows" about the 
higher dimensional origin of the underlying theory (the six dimensional (0,2) theory). 

lr This expression differs by a factor of 7r 3 from the original version of this paper. We thank C.G. 
Callan, A. Guijosa, K.G. Savvidy and 0. Tafjord for pointing out the wrong numerical coefficient. 
Also cqu. (3.12) had to be corrected, a factor of ^_ was missing. 
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3.3 't Hooft line 



In this section we calculate the energy between a monopole and an anti-monopole in 
the non-supersymmetric four dimensional theory obtained from the six dimensional 
(0,2) upon compactification. A related discussion for supersymmetric theories can be 
found in [41]. The supergravity background is the same as in the previous section. 
Namely, it is given by (3.8). As was explained in the previous section at large 
distances the effective theory is four dimensional along xi, x 2 , £3, X4. This theory 
being related to YM should contain monopoles. The string theory realization of the 
monopole is of a D2-brane ending on the D4-brane. The D2-brane is wrapped along 
x so from the point of view of the four dimensional theory it is a point like object. 
The action of a D2-branes is 

S = — / drda 1 da 2 e-*V^th, (3.21) 

where h is the induced metric. 

For our D2-brane which is infinite along one direction and winds the Xq direction 

we get 

Y 



9ym J \ -^1 
where we have used (3.9). Note the 1 / 9ym factor which is expected for a monopole. 
The distance between the monopole and the anti-monopole is 

U l J 2 h v /( y 3_i )(y3 _ 1 + e) 

where e = 1 — (Ut/U ) 3 . The energy (after subtracting the energy corresponding to 
a free monopole and anti-monopole) is 

2U 



dx< 



TT3 _ 7T3 

d x U* + U —^, (3.22) 



L = 2- 



j dV (3-23) 



E 



(2tt)3/V 



YM 



dy 



2(U T - U ) 

+ (27r) 3/2 • V-^) 



We would like to study the system in the region where LT 3> 1. At that region the 
energy turns positive which means that it is energetically favorable for the system to 
be in a configuration of two parallel D2-branes ending on the horizon and wrapping 
x which corresponds to zero energy (after the subtraction). So in the "YM region" 
we find no force between the monopole and the anti-monopole. Which means that 
there is a screening of the magnetic charge. 



4. Entropy and (No)-phase transitions in theories with max- 
imal supersymmetry 

In the conformal cases a phase transition cannot take place at finite temperature. 
The reason is that in conformal theories there is no scale and, therefore, a finite 
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T c cannot appear in the theory. Put differently, when compactifying the Euclidean 
time direction the radius of compactification cannot be considered as large or small 
simply because there is no scale in the theory to compare it with. In [46] SYM on 
S* 3 x S* 1 (rather than i? 3 x S* 1 for a superconformal theory in 3 + 1 dimensions) was 
studied. The radius of the S* 3 serves as a scale in the theory. Thus a phase transition 
is possible for this theory. Indeed, it was shown in [46] that a phase transition occurs 
at finite temperature. The phase transition was manifested through a change in the 
sign of the specific heat of the AdS black hole [43]. 

In this section we would like to consider the non-conformal theories with maximal 
supersymmetries which were studied from the supergravity point of view in [10]. 
These theories, being non-conformal, contain a scale which depends on gyn and N. 
Therefore, a phase transition at finite temperature in these theories might occur. At 
first sight it seems that phase transitions should appear in a natural way for the 
non-conformal theories. The reason is that different descriptions of the theories are 
valid in different energy regions [10] . However, analyzing the entropy of the systems 
describing the theory in these different regions, we do not detect any discontinuity 
in the entropy and the specific heat. This is an evidence that no first order phase 
transition occurs in these theories. Two comments are in order: (i) We cannot 
exclude possible jumps in the numerical prefactor of the entropy which, if it exists, 
will show a phase transition, (ii) We cannot exclude higher order phase transitions 
between regions other than the SYM and supergravity regions. Such transitions are 
known to occur in various cases as we enter into the eleven dimensional M-theory 
regime and are associated with "localization". An example of this phenomenon is 
the 2+1 dimensional theory [10]. 

As was reviewed in section 1 there are energy regions which are common to all 
p. These are the SYM region and the ten dimensional supergravity region. There 
are other regions whose description depends on p [10]. We shall first show that first 
order phase transitions do not occur in the common regions. Then we shall consider 
the other regions. 

4.1 SYM <-> lOd supergravity 

In this section we consider the SYM region, the lOd supergravity region and the 
transition between the two regions. Similar results using a slightly different language 
have been found in [48]. 

Perturbation theory in SYM is valid as long as g e ff « 1. In this regime the 
interactions between the gluons can be neglected and the free gas approximation can 
be used. For a free gas in d + 1 dimensions we have 

E~nVT d+ \ S~nVT d , (4.1) 

where n is the number of massless fields. So in our case n ~ iV 2 and hence we get 

Sym ~ Ar2/(P+iVV(P+i)£P/(P+i). ( 4 . 2 ) 
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The supergravity solution can be trusted as long as l/g e ff *C 1. (There is also a 
lower bound which we shall introduce in the next subsections) . The thermodynamics 
of the system in this region is defined by the area of the horizon of the supergravity 
solution which is [10], 

S sg ~ g^yfNEW^VW^. (4.3) 

First we would like to see whether the two entropies (4.2), (4.3) are of the same 
order at the transition region, g e ff ~ 1. Since g 2 ^ = g YM NU p ~ 3 the transition is 
at U ~ (gYM N Y /{3 ~ p) ■ The relation between U and E, V is [10] U J - p ~ gy M E/V. 
Therefore, the transition should take place at 

E = VN— g Y £\ (4.4) 
Indeed one can check that precisely at that point one gets 

S YM ~ S sg ~ l/^^pTV^i. (4.5) 
Let us consider now the specific heat c = ||? . From (4.2) and (4.3) we get 

Cy M ~ S'rAf, Csg ~ ^sg, (4.6) 

and hence at the transition cym ~ c sg . Moreover, it is clear from (4.6) that the 
specific heat is positive for any temperature and hence there is no first order phase 
transition in the SYM and/or lOd supergravity regions for any p. 

4.2 Dl-brane 

Eq.(4.3) with p — 1 cannot be trusted all the way to the IR limit (U — > 0). For 
U < gYM the proper description is via orbifold conformal field theory [47]. Before 
we discuss that region we should note that for U ~ gyuN 1 ^ the correct description 
is by the S-dual system. Since the entropy is defined in the Einstein frame and the 
Einstein metric is invariant under S-duality the S-dual description yields the same 
thermodynamics. 

In the region U < gYM the entropy should be calculated in terms of an orbifold 
conformal field theory. Like in the SYM case we consider the entropy in the free 
theory limit. The expression we find is a good approximation up to T ~ gym/VN 
[10] which corresponds to U ~ gYM- 

The maximal entropy is obtained when the configuration is that of one long 
string whose length is L = NL where L is the size of the system [45]. Therefore, we 
have 

E ~ LT 2 , S ~ LT (4.7) 

which gives 

S orb ~ VNLE. (4.8) 
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Since the transition occurs at U ~ qym and since U 6 ~ g YM EjL we find that at the 
transition region they are of the same order, 

Sorb ~ S sg ~g YM VNL. (4.9) 

Note that c orb ~ S or b and there is no first order phase transition at any finite tem- 
perature for SYM in 1 + 1 dimensions. 

4.3 D2-brane 

Again, like in the previous section, eq.(4.3) with p = 2 cannot be trusted all the way 
to the IR limit. At the point where [10] U ~ 9ym the correct description becomes 
the conformal theory of coinciding M2 branes with SO (8) R-symmetry. Therefore, in 
the region U < Qym the entropy is due to the area of a collection of N near-extremal 
M2-branes [44] 

S M2 = VNV 1/3 E 2 / 3 . (4.10) 
One can easily check [10] that near the transition region U ~ g YM 

Sm2 ~ S sg ~VNVg YM . (4.11) 

Note that cm2 ~ Sm2 and hence there is no first order phase transition that can be 
detected by the entropy . However, it was pointed out in [10] that there is a phase 
transition between a translation (along xu) invariant solution and a localized one. 
This type of phase transition, very likely to be a higher order one, does not show 
itself in the entropy or specific heat which are continuous. 

4.4 D4-branes 

Eq.(4.3) with p = 4 cannot be trusted all the way to the UV limit, U — > oo. The 
correct description for U > iV 1 / 3 /g YM is via M5-branes wrapped along x 10 [10]. The 
entropy is [44] 

S M5 ~ VNV 5 1/6 E 5 / 6 (4.12) 

where V 5 = V^ttRiq is the volume of the M5-brane. Since R w = gy M /(27r) 2 (4.12) 
is the same as (4.3) for p — 4. We conclude that there is no phase transition for 
the (0,2) theory on a circle. Note that for the six dimensional (0,2) theory a phase 
transition at finite temperature cannot occur simply because the theory is conformal 
and hence there is no scale. What we show above is that there is no first order phase 
transition even if one introduces a scale via compactification. 

4.5 D5-branes 

In the UV there is a transition from the D5-brane solution to the NS-fivebrane 
solution [10]. The transition is via S-duality. Since S-duality does not change the 
Einstein metric and the entropy is one quarter of the area in the Einstein frame the 
entropy is intact. 
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4.6 D6-branes 



The correct description of D6-branes in the UV region, U > -573-, is via Schwarzschild 

9ym 

black hole sitting on the A N _i singularity [10]. To calculate the entropy one should 
calculate the area of the horizon taking into account the identifications (which 
add a factor of 1/N to the usual Schwarzschild result) 



NE3 ' 2 *™, (4.13) 



V 

This coincides with eq.(4.3) for p = 6 and hence there is no first order phase transi- 
tion. 



5. Wilson loops and (no) phase transitions 

In this section we compute the space-time Wilson loops for the field theories with 
maximal supersymmetries both at zero and finite temperature. In these calculations 
we see no trace of phase transitions in the supergravity region. This is in accordance 
with what was derived in the previous section. 

In the four dimensional case, considered in [39], it was found that there are 
two regions. For TL < 1 we found a Coulomb behavior while for TL ^> 1 the 
configuration of two parallel strings which end on the brane is energetically favored 
and, therefore, the force between the quark anti quark vanishes. What we find for 
the non-conformal cases is exactly the same behavior. This is surprising because 
the non-conformal theories contain scales associated with the dimensionful qym in 
the relevant dimension. This supports, therefore, the conclusion of the previous 
section that there is no phase transition in field theories with maximal number of 
supersymmetries in the supergravity region. 

Before we calculate the Wilson loop with temperature let us first start with the 
Wilson line at zero temperature. 

5.1 Wilson lines at zero temperature 

In [37, 38] a stringy prescription for the computation of the expectation value of the 
spatial Wilson line of M = 4 was proposed. This expectation value determines the 
dependence of the energy E of a quark anti-quark pair on the distance between the 
quarks L. The goal of this section is to extract this dependence for non-conformal 
theories with 16 supercharges. In fact, one such case, the Af = 8 U(N) SYM in 
2+1 dimensions was analyzed in [37]. It was found that E ~ ( ^ Y ^2 — )^^- It was 
further shown [37] that this result, which is valid in the supergravity regime, is glued 
smoothly with the perturbative SYM result, which is valid in the UV. In the very 
IR there is a third regime where the correct description is in terms of the conformal 
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theory on coinciding M2 branes, but we will not attempt to treat this here. We 
address now the generalization of this result to other Dp-branes systems. 

The worldsheet action that corresponds to the metric G M n of Dp-branes takes 
the form 



S = — - — - I drdaJdetG MN d a X M d a X M 
2-Ka J v 

= L.J dxy/(d x Uy + IP-p/R^ . (5.1) 



where R p = (g YM d p Ny 

By repeating the procedure of [37] one finds 

L = 2^)^ f -^L_ i ~ 9YM VNUt^ (5.2) 

The energy of a system of a quark anti-quark is (for p ^ 5) 

E S9 ~ -(^)V(5-ri (5.3) 

the case of p = 5 is described separately below. This equation rests on the validity 
of the lOd supergravity description and hence it is valid as long as the range of U we 
integrate over is in the supergravity region defined in section (2). This has a double 
implication: (i) The minimal value of U, Uq, has to be greater than the lower bound 
stated in (2.4), (ii) since one cannot integrate up to U — oo, to ensure a reasonable 
approximation one has to demand that x(U u b) — L/2 <^ L, where U u b is the upper 
bound of U. For instance for p < 3 the consequence is that (5.3) can be trusted only 
for 

{g 2 YM Nyi^N^-W-rt « U « (g 2 YM Ny/(*-r) . (5.4) 

In certain cases discussed below we show that in fact the range can be extended. 

Next we show, as in the 2 + 1 case [37], that this lOd supergravity result is of 
the same order as the perturbative SYM result at the transition region. 

Eym ~ g YM NL 2 - p . (5.5) 

The transition occurs at U ~ (g YM N) 2 ^ 3 ~ p \ Since both energies (eq.(5.5) and 
eq.(5.3)) and the transition region depend only on g YM N, with no separate depen- 
dence on iV and g Y M, it is guaranteed on dimensional grounds, that at the transition, 
(5.3) 

E sg ~ E YM ~ 9 YM N^ ■ (5.6) 

To describe the behavior of the various systems in other regions of energy we need 
to do it for each p separately. 
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5.2 Dl-brane 



At first sight, according to (2.4), it seems that (5.3) with p — 1 is valid only in the 
region gyAiN 1 ^ <C U <C gyM^fN . In fact the domain of validity of (5.3) can be 
extended beyond the lower limit Uq < QymN 1 ^ where the dilaton becomes large [10]. 
In the region gyu <C Uq <C QymN 1 ^ 6 the proper description is via the S-dual system. 
Namely, the role of the "quark-anti-quark" is played by a Dl-string (instead of a 
fundamental string) and the background is of a collection of N fundamental strings 
(and not N Dl-branes). The action of a Dl-string in the background of fundamental 
strings is the same as the action of Fl-string in the background of Dl-branes. The 
reason is that due to the difference between their tension the Dl-string action is that 
of the fundamental string action multiplied with a factor e _9 \ On the other hand 
the S-dual metric contains a factor of so that these factors are canceled and one is 
left with the same expression for the energy of the system. Only at very low energies 
U < Qym [10] the orbifold conformal theory [47] takes over. 

We conclude, therefore, from the supergravity description, that for SYM in 1 + 1 
dimensions in the region 

IN , N 

, < L < , (5.7) 

^9ymN \j QymN 

the energy between a "quark anti-quark" pair is 



E ~(f%^ 1/4 - (5 ' 8) 

5.3 D4-brane 

Beyond the region where the ten dimensional supergravity is valid, for U ^> -t-^tf, 

9 Y m" 

the system is described by M5-branes wrapped on i? 10 [10]. In the case of M5 branes, 
which was addressed in [37], the role of the "quarks" is played by M2-branes wrapped 
on R w . The expression for the energy deduced from the "Wilson surface" in the six 
dimensional (2,0) theory is [37] 

where L is the length of the boundary of the M2-brane on the M5-branes. In our 
case since the M2-brane wrapped the x±o direction L' = 2-irRio ~ Qym- Thus we 
obtain the same result as the lOd supergravity expression (5.3) yields for p — 4. 
Thus, in the UV region L < Ng YM the energy is 

E ^- 9 1mE. (5.10) 
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5.4 D5-brane 

In the five-brane case there are only two regions. The SYM region and the 10D 
supergravity region. Like in the Dl-brane case the lOd supergravity region has to be 
divided into two region. But, again like in the Dl-brane case, the two regions, being 
related by S-duality, yield the same result. The solution of U(x) takes for p = 5 the 
following form U(x) = Uo/cos(x/R 5 ) which implies that L does not depend on Uq\ 

L = (2ir)- 3/2 g YM VN ~ VrtN. (5.11) 

This can also be seen also from (5.2) with p — 5. 

Let us consider the emerging physical picture. Suppose we start at large distances 
where we can trust the SYM calculation. The energy of the quark anti-quark is 
E ~ —g\ M N jL z . For L < gyM\fN we find that g e ff becomes larger then 1 and the 
supergravity description takes over. Then L is fixed to be gyM^/N and one cannot 
decrease L further. The fact that for L < gyAiVN there is no classical geodesic 
should be interpreted according to [46] as having a zero value of the Wilson loop 
which implies an infinite potential. However, a non-trivial dependence of U on L, 
and hence a non-zero Wilson loop, may emerge from a semi-classical quantization of 
the system with certain collective coordinates. 2 At this point we cannot resist to 
speculate that perhaps this "classical" minimal distance is related to the existance 
of a non-locality scale [42] . 



6. Wilson lines at finite temperature 



Next we address the large N behavior of the Wilson line of the non-conformal theories 
described in the previous section at finite temperature. As in the discussion of the 
p = 3 case this corresponds to a description in terms of a non-extremal supergravity 
Dp-brane [10] which translates into the following worldsheet action 

T 



S = — / dxy/(d x U) 2 + {W-p - U 7 T - p )/W-v 



(6.1) 



The determination of the E as a function of L,T, g YM y d p N follows the same steps 
as those leading to (5.3) subjected to the bounds on L mentioned above. 
The expression for the energy is 



E = 



Un 



dy 



' y/y?-r-l + e \ 
y/y<-P - 1 J 

where e = 1 — (Ut/Uq) 7 ~ p . For the length we find 
L = 2R£-v» 2 ut 5)/2 V~ej i 



+ 



dy 



(U T - U ) 



TT 



J(y7-P-l)(y7-P-l +e ) 



(6.2) 



(6.3) 



2 We thank J. Maldacena for pointing this to us. 
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At small temperatures where Ut <C U one finds the following result 



E (T,L) = co( 9 -^y^ 



1 + c(T _£_ ) ^)/(^) ) 



(6.4) 



The computation of the Wilson line for arbitrary temperature was performed for 
the case of p = 3 in [39]. The results for 1 < p < 4 are similar to that found in 
[39], for TL <C 1 the quark anti-quark pair is connected by a string and the energy 
goes as E ~ — £ 2 /(p-s) just like in the zero temperature case. For large TL > 1 a 
configuration of two parallel strings going from U = Ut to U = oo is energetically 
favored and the force between the quarks vanishes. The length where the string 
breaks is c/T where c is a dimensionless constant of order 1 (which does not depend 
on N, like in the 4d case [39]). We conclude that also in the non-conformal cases we 
find no phase transitions in the supergravity regime at finite temperature although 
there is a dimensionful parameter in the theory (in contrast to the conformal p=3 
case [39]). Maybe one has to choose a compact space to find a phase transition 
in field theories with maximal supersymmetries. For p=3 this has been done in 
[46] where it was found that the four dimensional theory on S* 1 x R 3 has no phase 
transition whereas the theory S 1 x S 3 exhibits a phase transition from confinement 
at small temperatures to deconfinement at high temperatures (the S 1 denotes the 
compactified Euclidean time direction with period (3 = T _1 ). But since in the 
non-conformal case the space-time is not a direct product of AdS p+ i times some 
compact space and we do not know the metrics that would correspond to spaces 
with boundaries of the form S* 1 x S p we have to leave the study of compact field 
theories to the future. 

The discussion in this section was restricted to the supergravity and Yang-Mills 
regimes but one may also try to connect to other regimes [10] e.g. in the p=2 case 
one needs the M2 brane description in the IR. 
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